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Abstract 

Predicting turbulent flame speed (the large time front speed) is a 
fundamental problem in turbulent combustion theory. Several models 
have been proposed to study the turbulent flame speed, such as the 
G-equations, the F-equations (Majda-Souganidis model) and reaction- 
diffusion-advection (RDA) equations. G-cquations describe the front 
motion law in the form of local normal velocity equal to a constant 
(laminar speed) plus the normal projection of fluid velocity. In level set 
formulation, G-equations are Hamilton- Jacobi equations with convex 
(L 1 type) but non-coercive Hamiltonians. In the first part of this 
paper, we show that flow induced strain reduces front speeds of G- 
cquations in periodic compressible and shear flows. The F-equations 
arise in asymptotic analysis of reaction-diffusion-advection equations 
and are quadratically nonlinear analogues of the G-equations. In the 
second part of the paper, we compare asymptotic growth rates of the 
turbulent flame speeds from the G-equations, the F-equations and the 
RDA equations in the large amplitude (A) regime of spatially periodic 
flows. The F and G equations share the same asymptotic front speed 
growth rate; in particular, the same sublinear growth law lpg ^ holds 
in cellular flows. Moreover, in two space dimensions, if one of these 
three models (G-cquation, F-equation and the RDA equation) predicts 
the bending effect (sublinear growth in the large flow), so will the 
other two. The nonoccurrence of speed bending is characterized by 
the existence of periodic orbits on the torus and the property of their 
rotation vectors in the advective flow fields. The cat's eye flow is 
discussed as a typical example of directional dependence of the front 
speed bending. The large time front speeds of the viscous F-equation 
have the same growth rate as those of the inviscid F and G-equations 
in two dimensional periodic incompressible flows. 
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1 Introduction 



Turbulent combustion is a complex nonlinear and multiscale phenomenon 
|30j . A comprehensive physical-chemical modeling requires a system of 
reaction-diffusion-advection (RDA) equations coupled with the Navier-Stokes 
equations. For theoretical understanding and efficient modeling of the tur- 
bulent flame propagation, various simplified or phenomenological models 
have been proposed and studied. Most notably, these models are passive 
scalar reaction-diffusion-advection equations (RDA) and Hamilton-Jacobi 
equations (HJ), as documented in books [351 133 ES] and research papers 
[B El El [91 dH E21 EH EH EH E31 EH EH] to name a few. For simplicity of 
presentation, throughout this paper, we shall consider the fluid velocity V 
as time-independent V = V(x). 

• (Scalar RDA model). The passive scalar reaction-diffusion advection 
equation for the temperature field is: 

T t + V(x) ■ DT = dAT + — f(T), 

T(x,0) = T (x), x E R n , (1.1) 

where T represents the reactant temperature, D is the spatial gradient op- 
erator, V(x) is a prescribed fluid velocity, / is a nonlinear reaction function; 
d is the molecular diffusion constant, r r > is reaction time scale. The 
flow field V is known or statistically known. For an isothermal reaction, the 
scalar is a reactant concentration however we shall still use T to denote it. 
The common form of the reaction function is f(T) = T(l — T), so called 
Kolmogorov-Petrovsky-Piskunov-Fisher (KPP-Fisher); f(T) = T m (l — T) 
(m > 2, higher order KPP-Fisher); f(T) = eT E l T {\ - T) (E > 0), Ar- 
rhenius combustion nonlinearity; F(T) = 0, T G [0,^] U {1}, f(T) > 0, 
T G (9, 1), ignition combustion nonlinearity. KPP or generalized KPP comes 
from isothermal autocatalytic reaction-diffusion system with equal diffusion 



constants (or unit Lewis number), [SIES]. Equation (1.1) is well-known to 
admit propagating front solutions if the advection is absent (V = 0), [2]. 
This is the case when a car engine filled with liquid or gaseous fuel is ignited 
to generate a spreading flame. Turbulent combustion concerns with the set- 
ting of flame propagation when the fuel is stirred on a broad range of scales 
for the purpose of speed enhancement and waste gas reduction. Though 
the flame front will be wrinkled by the fluid velocity, its average location 
eventually moves at a steady speed st in each specified direction, the so 
called "turbulent flame speed" . The prediction of the turbulent flame speed 
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is a fundamental problem in turbulent combustion theory |351 [3"T| [30] . For 
KPP nonlinearity, it is known HI [2T)l [3"5] that st is given by a variational 
principle on the large time growth rate of a viscous quadratically nonlin- 
ear Hamilton- Jacobi equation (QHJ). More precisely, considering compactly 
supported initial data T(x,0), then for each direction e and wave number 
A > 0, let £ e (A) be the principal Lyapunov exponent of the linear advection- 
diffusion equation: 

4>t = d A(f) + (-2 dXe- V(x, t)) • D<j> + [d A 2 + A e • V (s, y) + r" 1 / 7 ^)] 0, 

with initial data <^(x, 0) = 1. The function it = In 4> then satisfies the viscous 
QHJ: 

u t = dAu + d\Du\ 2 + (-2d\e-V(x,t))-Du + d\ 2 + \e-V(x,y)+T- 1 f'(0). 

Under suitable stationarity and ergodicity condition of the flow field [26] . 
the following limit exists almost surely and is independent of x: 

£ e (A) = lim -u(x, t). 

t— >+oo t 

The turbulent front speed along in e is: 

a deterministic quantity at large scale. In particular, if V is a periodic flow 
field, then 

£ e (A) =r- 1 f'(0) + H*(Xe). 
Here for p S K n , H*(p) is given by the following cell problem 

-dAw + d\p + Dw\ 2 + V(x) ■ (p + Dw) = H*(p) 

for w G C°°(T ?1 ), where T n is the n-dimensional flat torus. 

• (Majda-Sougandis model). Turbulent combustion might involve small 
turbulent scales which are no less than the reaction time scale in the thin 
flame model. In [21] . the velocity field V is space-time periodic and in a scale- 
separation form V = V(x,t,e~ a x,e~ a t), d = ed, r r = er r , a G (0,1]. For 
simplicity of analyzing upscaling, we consider V as time-independent with 
no integral scales, i.e V = V(e~ a x). The limiting behavior of T = T e is 
[21] : lim e ^o^ e = locally uniformly in {(x,t) : Z < 0} and T e — > 1 locally 
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uniformly in the interior of {(x, t) : Z = 0}, where Z G C(M n x [0,+oo)) is 
the unique viscosity solution of the variational inequality 



max(Zi -H(DZ) -T~ l f{0),Z) = 0, (x, t) x K n x (0, +oo), (1.3) 

with initial data Z(x,0) = in the support of T(x,0), and Z(x,0) = — oo 
otherwise. The set T t = d{x G W 1 : Z(x, t) < 0} can be regarded as a front 
which moves with normal velocity 

v n = F{n) 

where n denoting the normal vector pointing to the propagation direction 
and 

A>0 A 

The effective Hamiltonian H = H(p) is defined as a solution of the following 
cell problem: for each p G M. n , there are a unique number H(p) and a 
function F(x) G C - 1 ^™) such that 

- a{a)dAF + d\p + L>F| 2 + V(x) • (p + DF) = (1.4) 

where a(a) = if a G (0, 1), a(a = 1) = 1. The turbulent flame speed along 
the unit direction e is: 

s T {e)=F{e). (1.5) 



Note that when a = 1, (1.5) coincides with rtl.2h. Hereafter, we shall refer 



to (1.4) as F-equation if d = 0, and viscous F-equation if <i > 0. 



• G-equation (Level Set) model. The RDAs (1.1) and in particular the 
KPP are first principle equations, yet they are limited to the unit Lewis 
number regime. Another approach in turbulent flame modeling is the level 
set formulation of interface motion laws in the thin interface regime. The 
simplest motion law is that the normal velocity of the interface (y n ) is equal 
to a constant s/ (the laminar speed) plus the projection of fluid velocity along 
the normal n. See Figure 1. The laminar speed is the flame speed when 
fluid is at rest. Let the flame front be the zero level set of a function G(x, t), 
the burnt region is G(x, t) < 0, and the unburnt region is G(x, t) > 0. The 
normal direction pointing from the burnt region to the unburnt region is 
DG/\DG\, the normal velocity is —Gt/\DG\. The motion law becomes the 
so called G-equation, a popular model in turbulent combustion [351 GO] : 

G t + V(x)-DG + Sl \DG\=0. (1.6) 
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Figure 1: Illustration of G-equation (level set) model. 



Chemical kinetics and Lewis number effects are all included in the laminar 
speed si which is provided by a user. Formally, under the G-equation model, 
for a specified unit direction p, 

s T (p) = - lim (1.7) 

T-»+oo I 



HereG(x,t) is the solution of equation (1.6) with initial data G(x, 0) = p-x. 
We remark that in combustion literature, the level set function is defined as 
G(x, t) > in the burnt region and G(x, t) < in the unburnt region (pp 91- 
92 of [30J). The resulting G-equation is G t + V{x, t)-DG = s t \DG\. We shall 
work with the form in ( |1.6[ ) instead. The existence of st has been rigourously 
established in |37| and [6] independently for incompressible periodic flows. 
And st is the effective Hamiltonian of the following cell problem 

si\p + DG\ + V(x) ■ (p + DG) = H{p) = s T . (1.8) 

Here H(p), the effective Hamiltonian, is the unique number such that the 
above equation admits periodic approximate solutions. See the Appendix 
for some basic facts about viscosity solutions and cell problems. The for- 



mal analysis of (1.6) and st is also performed in the framework of renor- 
malization group methods [331 EH EH] . See also [29J on a spectral closure 
approximation, and [T7] for a numerical study of G-equation in comparison 
with a combustion system modeling thermal-diffusive instabilities of free- 
propagating premixed lean hydrogen-air flames. Though G-equation is a 
phenomenological model, it is more flexible in that many factors influencing 
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front motion can be incorporated into s/. For example, the strain effect of a 
turbulent fluid flow is modeled by extending si to si + d'rt ■ DV ■ n, where 
d is the Markstein length which is proportional to the flame thickness. The 
G-equation with flow induced strain is [30] : 

DC 

G t + si\DG\ + V(x, t)-DG + d—- -DV-DG = 0. (1.9) 

\L>G\ 



Then formally, st is also given by (1.7), where G(x,t) is the solution of 



equation (1.9) with initial data G(x, 0) = p ■ x. So far we are not able 
to prove the existence of st except for some simple situations like the one 
space dimensional (Id) compressible flow and the incompressible shear flow. 
It is conjectured by some experts [31J in combustion theory that the strain 
term will slow down flame propagation. Theoretically it is hard to verify 
this for general flows. In section 2 on this paper, we confirm this conjecture 
for Id compressible flow and the shear flow. Besides the strain effect, the 
flame stretching mechanism also includes the curvature effect. One such 
model proposed in [30] is to replace si by si(l — dn). Then the G-equation 
becomes 

G t - ds l K\DG\ + si\DG\ + V(x, t) ■ DG = 0. (1.10) 
Here d is the Markstein length and k is the mean curvature of the flame 



front, i.e, k = div (tRt^t)- The curvature G-equation (1.10) is very difficult to 



■ \DG\ 

analyze. To obtain some ideas of the diffusion effect, a na" 



;ural simplification 



is to change the mean curvature term k to AG. This leads to the viscous 
G-equation 

G t - dsiAG + si\DG\ + V{x, t) ■ DG = 0. (1.11) 

The above viscous G-equation also serves as a basic model to understand 
the numerical diffusion effect introduced in the numerical computation of 
equation (|1.6[). For the viscous case, st = ST(p,d) is given by the large 



time limit ( |1.7[ ), where G(x,t) is the solution of equation (1.11) with initial 
data G(x, 0) = p ■ x. It is also the effective Hamiltonian of the following cell 
problem 

-dsiAG + s t \p + DG\ + V{x) ■ (p + DG) = s T {p, d) = H(p, d). 

Here H(p, d) is the unique number such that the above equation admits 
periodic solutions. The most general G-equation is to combine both the 



strain effect and the curvature effect, that is si in the basic G-equation (1.6) 
is replaced by si(l — dn) + dlt ■ DV ■ it. 
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Though various passive scalar models as shown above have been pro- 
posed to study st, their predictions may be potentially different or some- 
times asymptotically identical. It requires delicate analysis to understand 
these subtleties. The goal of this paper is to analyze and compare qualita- 



tive and quantitative properties of st in G-equations (1.6)-(1.9), F-equation 



(1.4)-(1.5) and the RDA equation (1.1 )-( 1.2 ) for steady (time-independent) 
periodic flows (V = V(x)) with mean equal to zero. The s^'s are compared 
in terms of different nonlinearities, and the flow induced strains. 

Outline of the paper. To simplify notations, throughout this paper, we 
normalize s; = d = r r = 1. 



Section 2: We analyze G-equation (1.6)-(1.9) with and without the strain 



effect. We show that the strain term slows down the propagation speed st 
in one- dimensional compressible flows and in shear flows. 



Section 3: The comparison of the F-equation (1.4)-(1.5) with the G-equation 



(1.6) had been studied in [IT] for periodic shear flows. Their computations 



showed that G-equation model always predicts slower turbulent flame speeds 
than the F-equation model for the shear flow. We will prove this result for 
general flows as long as /'(0) is no less than a threshold value. We also 
investigate the asymptotic behavior of st from the reaction-diffusion equa- 



tion (1.1)-(1.2), the G-equation (1.6) and the F-equation (1.4)-(1.5) when 



V(x) is scaled to AV(x), i> 1. We prove that the asymptotic growth rate 
lim^^+oo ^ is the same for G and F equations. The limit is given by an 
inf-max formula which is also equivalent to a variational formula involving 
the invariant measures of the flow x = V(x) based on the weak KAM theory. 
We also provide a necessary and sufficient condition for the bending effect 
(sublinear growth in A) in two space dimensions (n = 2) in terms of the 
periodic orbits of the dynamical system x = V(x) on the two dimensional 
torus and their rotation vectors. We show that when n = 2 the asymp- 
totic growth rate of st from the RDA equation is qualitatively the same 
as st from the G and F equations. In particular, these three models will 
predict the bending effect simultaneously. See Theorem 3.3 for the precise 
statement. As an application, the growth laws of st as a — > +oo predicted 
by these three models are completely identified for the cat's eye flow and 
the cellular flow. In particular, for cellular flows, the i og ^) g rowth law of 
st in G-equation (1.6) is also verified for the F-equation (1.4)-(1.5). These 



two flows appeared as representative flow examples in related dynamo and 
convection-enhanced diffusion problems [T3] . 

Section 4: We prove that in two space dimensions, the effective Hamilto- 
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nian from the viscous F-equation gives the same asymptotic growth rate as 



st in A predicted by the F-equation (1.4 )-( 1.5 ) and the (inviscid) G-equation 



(1.6). It follows that the effective Hamiltonian of the viscous F-equation 
grows sublinearly in large cellular flows. On the other hand, it has a lower 
bound of OiA 1 / 4 ), see [21] . Moreover, we also show that, in 2d, changing 
the order of sending the flame thickness or reaction time scale to zero and 
sending the turbulent intensity to infinity does not change the asymptotic 
growth rate of the front speed st- 

Appendix: For reader's convenience, we will review some basic facts about 
viscosity solutions of Hamilton-Jacobi equations and their cell problems. 

Acknowledgements: The work was partially supported by NSF grants 
DMS-0911277 (JX) and DMS-0901460 (YY). 

Assumptions and notations: Throughout this paper, we assume that 
nonlinear function T — > f(T) is of KPP type, i.e 

/eC^ao,!]), /(0) = /(l) = 
and / is non- increasing on (1 — e, 1) for some e > 0. Also 

< f(s) < s/'(0) for s€ (0,1). 

Note a prototypical example is 

f(T) = KT(l — T). 

(1) Z n = {(xi,X2, ...,x n )\ Xi G Z}. T n denotes the n-dimensional flat torus, 
i.e, T n = R n /Z n . 

(2) / G C r {T n ) if / G C r (R n ) and it is periodic, i.e, f(x + v) = f{x) for 
x G W n and v G Z n . 

(3) Df is the gradient of /. If g G M n -> M. 

(4) Throughout this paper, we assume the velocity field V G C°°(T n ) and 
has mean zero, i.e, J T „ V dx = 0. 

(5) For a curve £ : M — > W 1 , £ is its image under the natural projection 
W 1 — > T n . See Figure 2 for a projected periodic orbit on the torus. 

If we scale V to AV: 

(6) (xa(p) the turbulent flame speed predicted by the G-equation equation 



model (1.6)-(1.8). 



(7) 7a(p) the turbulent flame speed predicted by the F-equation equation 



model fll.4D-fll.5D for a G (0, 1). 

(8) c*(A) the turbulent flame speed predicted by the RDA equation model 

ob-diib. 



7 



Figure 2: Schematic of the projection M n — > T n , and a periodic orbit on T n . 

2 G-equations with and without Strain 



The goal of this section is to show that the strain term in (1.9) slows down 
flame propagation for Id compressible flow and the shear flow. This is differ- 
ent from the diffusion effect which actually enhances the flame propagation 
for Id compressible flow [19J. Recall that the turbulent flame speed pre- 
dicted by the G-equation is the effective Hamiltonian associated with the 



cell problem (1.8). See the Appendix for some basic facts about viscosity 



solutions and cell problems. 

2.1 One-dimensional G-equation and Compressible Flow 

It suffices to look at p = 1. We write the flow velocity V(x) = v(x). For 
one-dimensional G-equation without the strain term, the turbulent flame 
speed is the unique number c such that the following equation (p = 1) 

\l + G'\+v(x){l + G') = c 

admits approximate periodic viscosity solutions. Let us assume that J v(x) dx 
to make c > 0. It has been proved in section 3 of [19] that 

JO if {x\ v(x) = -1} ^ 
C ~ T+k?) dxf 1 > 0. 

Now let us consider the G-equation with the strain term. Using the same 
method as in section 3 of |19] , we can show that there exists a unique number 
c such that the following equation 

(l + v')\l + G'\ + v{x)(l + G') = c 

admits approximate periodic viscosity solutions. The c is the turbulent flame 
speed predicted by the G-equation under strain, and is given by 

[O if {x\ v'[x) +v(x) = -1} ^ 

-l 



. (Jo l+v'{x)+v(x) dx ) > °- 
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Theorem 2.1 The flame speed in the G-equation with strain is no faster 
than that of the G-equation without strain: 

c> c. 

If c > 0, then "=" holds if and only if v' = 0. 
We first prove the following lemma. 

Lemma 2.1 Suppose that w G C 1 (T 1 ) and w + w' > 0. Then w > and 



f- 

Jo « 



1 1 



— dx < I dx. 

w J w + w 

The "=" holds if and only if w is a constant. 

Proof: Assume that w(xq) = min^i w. Then 

w(xq) = w(x ) + w'(x Q ) > 0. 



So w > 0. Now let s = Then = w ^ w , ■ Since s — s' > 0, Cauchy's 
inequality implies that 



1 / f 1 s 2 ( f 1 ^ 2 

(s — s) dx / dx > ( / s (ix 



o 



Note that L(s — s')dx = Jq sdx. Accordingly, 

/ dx > s dx. 

Jo s-s' J 

The equality holds if and only there exists A > such that 

s — s' = As. 

Since s is periodic, this implies that s' = 0. So w' = 0. □ 



Proof of Theorem 2.1 Suppose that c > 0. Then 1 + v + v' > and the 

above theorem follows immediately from Lemma 2.1 by setting w = 1 + v. □ 
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2.2 Strain Effects in Shear Flows 

Suppose that V(x,y) = (v(y),0) for (x,y) G R 2 . For (m,n) G M 2 , denote 
A(m,n) as the unique number such that the G-equation for the shear flow 

y 7 m 2 + (n + u') 2 + mv(x) = A = A(m, n) 

has a periodic viscosity solution. We also write A(m, ra) as the unique number 
such that the G-equation for the shear flow with strain term 

/ — ^ -, ttk m(n + u')v' , . - - . . 

y m z + (ra + H + mv{x) = A = A(m, raj. 

Yrra 2 + (ra + u') 2 

It is clear that 

A, A > |m| + maxmv (2-1) 
ji 

To ensure that both of them are nonnegative, we further assume that 

1 

v(x) dx = 0. 

Note that A(m, ra) and A(rra, ra) are the turbulent flame speeds predicted by 
the G-equation without and with strain effects respectively for the shear 
flow. 

Theorem 2.2 The flame speed under the strain of shear flows is no faster 
than that in shear flows without the strain. 

A(m, n) > A(rra, n). 

If m 7^ and A(m,ra) > \m\ + maxp m», then "=" holds if and only if 
v' = 0. 

We first prove several lemmas. 

Lemma 2.2 Given m,c£l and a > \m\, there exists a unique p > such 
that 

I — 9 9 C P 

\fm z + p z H ^^^^= = a. 

Y rat 2 + p 2 

Proof: Denote /(p) = y 7 m2 + P 2 H — 7 ? „ • Clearly there exists at least 
one po > such that 

f(Po) = ci- 
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Note that if f(p) = a and p > 0, then p + c > and 

ft,) = ( f;t m2 t,f > »■ 

So po is unique. □ 



Lemma 2.3 Suppose that s(x) and k{x) are C°° periodic functions. Given 
m/0,n/0,c6l and a>\m\ +max T i k, there exists at most one periodic 
viscosity solution u satisfying 

l — k -, tttt c(n + u f )s(x) , 

Vm 2 + (n + u') 2 H + = a. 

ym 2 + (n + n') 2 

Proof: Without loss of generality, we assume that n > 0. Clearly, u G 
W 1 ' 00 ^ 1 ). Due to Lemma 2.2 to prove that it is unique and C , it suffices 
to show that 

n + -u'>0 a.einT 1 . (2.2) 

We claim that w(x) = nx + -u(x) is increasing. In fact, if not, then w must 
have a local minimum point since w(x + 1) — w(x) = n > 0. Suppose w 
attains local minimum at xq. Then by the definition of viscosity solutions, 
we have that 

\m\ + k(xo) > a. 

This is a contradiction. Hence our claim holds. Note that w' can not be 
zero. Therefore (2.2) is true. □ 

Assume that s, k are smooth periodic functions. Fix m and n. For c > 0, 
we denote h(c) as the unique number such that the following equation 

y/m 2 + (n + h') 2 + < n + u ')-< x ) +k(x) = h(c) (2.3) 
W m? + (n + u') 2 

admits viscosity solutions. The existence of h{c) follows from |18| . 

Theorem 2.3 Suppose that J Q s(x) = and s(x) is not identically 0. Also 
we assume that 

s(x) = if k(x) = max k. (2-4) 
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If m, then 

(i) h(c) is Lipschitz continuous and min^+oo) h= \m\ + max T i k; 

(ii) h'(c) < for a.e c E (0, oo); 

(Hi) If h(c) > \m\ + max T i k, then h'(c) < 0; 

(iv) There exists Co such that h(c) = \m\ + max T i k, for c > cq. 

Proof: (i) and (iv). We claim that 

\h(ci) - h(c 2 )\ <L\ci-c 2 \ (2.5) 

for L = max T i In fact, suppose that u\ and u 2 are viscosity solutions 
of (2.3) for h(ci) and h(c 2 ) respectively. We may assume that both of them 
are C 1 . Otherwise, we just apply the routine arguments using the method 
of "doubling the number of variables" . Suppose that 

ui(x ) - u 2 (x ) = max(ui - u 2 ). 

T 1 

Then 

u x {xq) = u 2 (x ) 

and 

h( Cl ) - h{c 2 ) = (ci - C2)g(XQ)(n + U l (X0 )) 
Jm 2 + (n + u'^xq)) 2 

Hence (2.5) holds. Now owing to (2.4), it is obvious that 

h > \m\ +maxfc. (2.6) 

T 1 

Since s is not constantly 0, there exists xq G (0, 1) such that r = s(xo) < 0. 
Choose e > such that [xo — e, xo + e] G (0, 1) and for x £ [xq — e, xo + e] 

— r 

six) < — . 
v j - 2 

Let (ft be a C°° periodic function such that n + (p' > and within the interval 
[0, 1] 

supp(n + cj)') C (x - e, x + e). 
Now choose cq > such that 

cqt > 2 max(n + 4>'). 
T i 

Then it is clear that for all c > cq 

( I — t / 77v5 c(n + (f)')s(x) 

max{-y/m z + [n + cp') z H -j^^^^=^^^^= + fc(x)) = m + max k[x). 

T 1 y/m 2 + (n + ^') 2 T1 
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Suppose that u(x, c) is viscosity solution of equation (2.3) for c > cq. Con- 
sidering the place where u — <j> attains minimum, we derive that for c > Co 

h(c) < \m\ +maxfc(x). 



Combining with (2.6), (i) and (iv) follows. 

(ii) and (iii): Suppose that h is differentiable at c. If h(c) = \m\ + 
max T i k(x), then it is clear that h'(c) = 0. So let us assume that h(c) > 
\m\ + max T i k{x). Since h is continuous, when c is close to c, h(c) > \m\ + 
max T i k(x). Owing to L emma 2.3, let u = u(x,c) be the unique smooth 



we assume that n > 0. 



solution of equation (2.3) subject to udx = 0. Without loss of generality, 



hen n + u' > 0. Denote 



f[p) = \fm 2 + 



Then equation (2.3) can be rewritten as 



f{n + u ) + cs(x)f'(n + u ) + k(x) = h(c). 
Taking derivative with respect to c on both sides, we have that 

h'(c) = f'(n + u')u c + cs{x)f"(n + u')u c + f(n + u')s(x). 
Note that 

m 2 (p(x) + cs(x)) + p 3 (x) 



f'{n + u ) + cs(x)f"(n + u ) 



for = n + u' . Since > and h(c) — \m\ > max T i k, owing to (2.3), 
we have that 

p(x) + cs(x) > 0. 
Hence f'(n + u') + cs(x)f"(n + v!) > 0. Therefore 



f's(x) 



+ 



h'(c) 



f' + f"cs(x) f + f"cs(xy 
Taking integration on both sides over [0, 1], we obtain 

1 



h'(c) 



f + f"cs{x) 



dx 



six 



1 + a(x)s(x) 



dx, 
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where a(x) = ^ir > 0. Since 1 + a(x)s(x) > 0, Jq s(x) 
identically 0, 

s(x) 



and s(x) is not 



So for c close to c 
□ 



o l + a(x)s(x) 
h'(c) < 



dx < 0. 



Proof of Theorem 2.2 Case 1: If m = 0, then it is easy to see that 
A = A = |n|. 

Case 2: If n = 0, then A = A = \m\ +max T i mv. In fact, let u be periodic 
viscosity solution of 



y/m 2 + (u') 2 + 



u'mv' 



+ mi; = A. 



a/771 2 + (u') 2 

Suppose that u(xq) = min T i u. Then 



A < \m\ + mv(xcj) < \m\ + rn.ax.mv. 

T 1 



Combining with (2.1), A = \m\ + max T i mv. The proof of the equality 
A = \m\ + max-p mv is similar. 



Case 3: If mn ^ 0, then it follows immediately from Lemma |2.3| by 
choosing s(x) = mv' and k(x) = mv. 



3 Comparison of Asymptotic Growth Rate of st 
from G-equation, F-equation and RDA 

The main purpose of this section is to find quantitative and qualitative 
similarities between asymptotic growth rate of the turbulent flame speeds 



predicted by G-equation, F-equation and RDA. Theorem 3.3 implies that 
in two dimensions, these three models predict the bending effect simulta- 
neously. We also provide a simple necessary and sufficient condition of the 
nonoccurrence of the bending effect. The growth laws of st as a — > +oo 
predicted by these three models are completely identified for the cat's eye 
flow and the cellular flow. This partially answers a question posed at the 
end of [H]. Recall that we normalize d = si = r r = 1. Throughout this 
section, we assume that V{x) is divergence free and has mean zero. When 



14 



n = 2, this is equivalent to saying that there exists a periodic smooth stream 
function H such that 

V(x) = (-H X2 ,H X1 ). 

3.1 A Rough Comparison of St from G and F-Equations 

For A > 0, we scale V to AV. For p 6 M n , we denote aA(p) as st from the 
inviscid G-equation which is the effective Hamiltonian of the following cell 
problem 

\p + DG\+AV(x)-(p + DG) = a A {p). (3.1) 

We also denote (3a{p) as the effective Hamiltonian from the following cell 
problem of invisicd F-equation 

\p + DF\ 2 + AV(x) ■ (p + DF) = P A (p)- (3.2) 



We recall that oia(p) is the unique number such that equation (3.1) admits 



approximate periodic viscosity solutions and (5a (p) the unique number such 



that (3.2) admits periodic viscosity solutions. Note that ( |3.l| might not 
have exact solutions due to the lack of coercivity. Both a a (p) and (3a(p) 
can be given by inf-max formulas. Precisely speaking, 

a A (p)= inf max{|p + ZM +AV(x) ■ (p + D6)} (3.3) 

(>eC 1 (T n ) T" 

and 

Pa(p)= inf m&x{\p + DS 2 + AV(x) ■ (p + D6)} (3.4) 
Hence a^(p), Pa {p) <0(A). Since t 2 > t - \, 

pA>a A -\- (3.5) 
We write 7a(p) as st from the Majda-Sougandis model along direction p for 

/'(0) + fli(Ap) 



a G (0,1). Then by (1.5) 



7a(p) = inf 



A>o A 



Lemma 3.1 Denote r = 2y / /'(0). T/ien 

ra^Cp) < 7 a(p) </3a(p) + /'(0). 
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Proof: The right inequality is obvious by choosing A = 1. Let us prove the 
left inequality. Since > t — J, according to the inf-max formulas, 

^(p)>rax(p)-/'(0). 

T 

Combining with the degree 1 homogeneity of au(p) with respect to the p 
variable, the above lemma holds. □ 



Remark 3.1 Note that if we choose f'(0) > \ as in fllf . then ja(p) — 
(Xa{p)- This generalizes the computations in Jii] / for shear flows, which 
shows that the G-equation model always predicts slower turbulent flame speeds 
than the F-equation model. However if /'(0) < i, this conclusion is not true 
in general. Here is a simple example based on the shear flow. Let n = 2, 
V(xi,X2) = (v(x2),0) and p = (1,0). Then an easy calculation shows that 

a A (p) = 1 + A max v 

T 1 

and 

Pa(^p) = A 2 + AXm&xv. 

T 1 

Then 7^ (p) = Amax T i v + 2y / /'(0) < o.a(p)- 



3.2 Identical Asymptotic Growth Rate for G and F Equa- 
tions 

Throughout this subsection, for convenience, we drop the dependence of 
a Aip)i Pa{p) and 7a(p) on p. The following theorem says that a a/ A, Pa/ A 
and ^a/A have the same asymptotic limit. 

Theorem 3.1 Given p^W 1 . Denote 

Cn = inf maxjy(x) • (p + D6)\. 

Then 

,. ol A /3 A 7a 

lim — — = lim — — = lim — — = c„. 

A^+oo A A-s>+oo A A^+oo A 

In particular, G-equation and F-equation models predict the bending effect 
simultaneously. 
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Proof: The proof is simple. Owing to the inf-max formula (3.3) 

OiA . 

Now fix e > and choose (p t G C 1 (T ri ) such that 

max{T/(x) • (p + D(p t )} < e + c p . 

Then 



Hence 



< — max \ D(f) e \ + c p + e. 



,. OLA 

lim — - = On. 



The proof for (5a is similar. The proof for 7^4 follows from Lemma 3.1 □ 



Now an interesting question is when c p > occurs, i.e, there is no bend- 
ing effect. To this end, we employ the weak KAM theory to introduce an 
equivalent formula of c p . 



Lemma 3.2 



c p = max / p ■ V(x) da, 



where A is the collection of all Borel probability measure on T n invariant 
under the flow £ = V(£). 

Proof: Assume that (f) 6 C 1 (T n ) and a £ A. We have that 

/ p-V(x)da= p-V(x) + D(j)-Vda <max{V(x) ■ (p + D<f>)}. 



Hence 



Cp > max / p ■ V(x) da. (3-6) 



Denote the Hamiltonian 



H A (p,x) = ^\p\ 2 + V{x) - p. 



Then the corresponding Lagrangian is 



A 2 
L A (q,x) = — \q- V(x)\ . 
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Write H A (p) = Then H A (p) is a convex function of p. Now fix p 

and choose Qa £ 9Ha(p)- Here 8Ha(p) is the set of subdifferentials of Ha 
at p, i.e, 

cW A (;p) = {Q G M n | #V) > ff(p) +Q-{p -p) for all p' G M n }. 

Let /i^ be a Mather measure associated with the Lagrangian La with the 
rotation vector Q. Then ha is a Borel probability Lagrangian flow invariant 
measure on the phase space W 1 x T n which minimizes the functional 

L(q, x) dfjL 

"xT" 

among all Borel probability Lagrangian flow invariant measures fi on the 
phase space R n x T n subject to 



/ 

Je 



qd/j-A = Qa- 

"xT" 



See [23] for the existence of Mather measures. Also, suppose that F is a 
viscosity solution of 

j\P + DF\ 2 + V(x) ■ (p + DF) = H A (p). 

Then F is differentiable on the projection of supp(/u) to T n and 

p + DF(x) = — (q— V(x)) for (q,x) G supp^). 



See [13]. Due to Lemma |3.3[ 

\q - V(x)\ < o(l) for (q,x) G supp(/u), (3.7) 
where lim^^+oo o(l) = 0. Since Ha(0) = 0, we have that 

P-Qa> H A (p) > c p . 
Upon a subsequence if necessary, we may assume that 

HA \i weakly in R"xT" 



Suppose that a is the projection of \i on T n . Owing to (3.7), 

= lim^^+oo p- Qa>c p . 
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and 

supp// C {(q,x)\ q = V(x)}. (3.8) 

Now we show that a is flow invariant. In fact, since ha is Euler-Lagrangian 
flow invariant, for any (ft £ C 1 (T n ), 



/ 



q ■ D(j) dfiA = 0. 

"xT" 



Sending A — > +00 and using (3.8), we derive that 

/ Dcj) ■ V da = 0. 

Therefore a G A and 



/ p ■ V{x) da > c p . 



Combining with (3.6), Lemma 3.2 holds. □ 



Lemma 3.3 Suppose that F is a viscosity solution of (3.2), then 

DF(x) 



lim sup 

A-^+oo xG m 



0, 



A 

where Ai is the set where F is differentiate. 

Proof: Throughout this proof, C denotes a constant depending only on V 
and \p\. Choose F such that L„ F dx = 0. Denote w{x) = P ' x + F . Then w 
satisfies that 

|-Du?| 2 + V(x) ■ w(x) = — 

Step 1: we claim that 

lim w = uniformly in T n . 

In fact, since < 0(A), we have that \Dw\ < C. Upon a subsequence if 
necessary, we assume that 

lim w = h uniformly in T". 

Then h is a mean zero periodic viscosity solution of 

\Dh\ 2 + V(x)- Dh = 0. 
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Taking integration on both side, we obtain that 

/ \Dh\ 2 dx = 0. 

Hence h = 0. Our claim holds. 
Step 2: we claim that 



sup \Dw(x)\ < C ( —j= + ^/max 



/ 


F 


max 




T n 


~A 



A 

Choose xo £ T ra such that w is differentiable at xq. Let : (— oo,0) be 
the backward characteristics with £(0) = xq. Then 



w(x ) - w{t{-t)) = + \ [° \i - V(0\ 2 ds. 



A 2 '4. 

Also, £ satisfies the Euler-Lagrange equation 

a = dv(o ■ i - a - v) ■ dv(o 

and the equality 



C - v(0 



(3.9) 



(3.10) 



(3.11) 



Accordingly, |f| < C. Recall that w = Choosing t = 1 in <\3.9\), 

we deduce that 



,2- -C. 



\t-V(0\'ds<-(l + max\F\). 



So there is a to £ [0, 1] such that 



C. 



|£Ho) " ^(?(-*o))r < tC 1 + max|F|). 



By equation (3.10), 
d 



ds 



\i(s)-v(a(s))\ 2 <c\i(s)-v(as))\' 



This implies that 

|£(0) - V(m)\ 2 < e ct °\a-t ) - V(a-t ))\ 2 - 



Combining with (3.11 ), our claim holds. Lemma 3.3 follows from Step 1 and 
Step 2. □ 
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Definition 3.1 A smooth curve £ : M — > R n is called an orbit if = 
V(£(t)). Moreover, £ is called periodic orbit and T > is called a period if 
£ is an orbit satisfying that £(T) — £(0) is an integer vector. Furthermore, 

is called the rotation vector of £. 

Theorem 3.2 Given peW 1 . 
(i) For any orbit £, 

limsup P ' < c p . (3.12) 

In particular, if £ is a periodic orbit with period T , then 

„ > M^H(o)) 

(ii) There exists an orbit £ such that 

(Hi) When n = 2, there exists a periodic orbit £ with period T such that 

MM)-, 

where T is the period. 

Proof: (i) In fact, for any <j> G C^T"), 

Hence 

lim sup T ^ +oc £^ = lim sup T _^ +00 ±jj P ■ V (£(*)) eft 
< max T n + D<f>) ■ V} . 
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So fl3l2| ) holds 



(ii) Choose do € A such that 

/ p ■ V(x) dao = m & x / P • V(x) da = c p . 

By the Birkhoff Ergodic Theorem for measure preserving flows, there exists 
a o"o measurable function t/j such that for ao a.e x 

'ib(x) = lim - — - 

and 

p-F(x)dcr = / ifi(x)dao, 



where £ x = V(6e) and ^(0) = x. By (3.12) 



Hence for ao a.e x 



V ; ( x ) < / p ■ V"(x) doo- 
Jt 2 



(hi) Now let us assume n = 2, according to Poincare recurrence theorem, 
a.e. in x are recurrent points. Note that if x is recurrent, £ x must be 
periodic. The reason is simple. Assume that V(x\,X2) = (— H X2 , H Xl ) for 
some smooth stream function H. Then H is constant along £ x . Hence £ x 
must be periodic. 

The following is an immediate corollary. 

Corollary 3.1 Assume n = 2 and p£l 2 . Then 

c p > 

z/ and only there exists a periodic orbit £ wii/i period T such that 

p • (£(T) - £(0)) > 0. 

Lemma 3.4 Suppose that £ : [0, T] I 2 is a periodic orbit and T is the 
period. Denote Q = ^ T ^~^ ^ , Assume that Q / 0. TTien 
(%) c p >0ifp-Q> 0; 

f«J a A (p), /3a(p), 7a(p) < 0(1) forallA>0ifp-Q = 0. 
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Proof: (i) follows immediately from Corollary |3.1| Let us prove (ii). Through- 
out this proof, for i £ R 2 , ^ denotes the orbit satisfying ^(0) = x. Owing 



to (3.5) and Lemma 3.1, it suffices to prove for /3a- The strategy is to 



construct a suitable function 5 S W 1,00 (M. 2 ) satisfying 

V ■ (p + DS) = 0. 

Denote x = £(0). Choose a stream function H such that V = V^H = 
(— H X2 , H X1 ). Since H is constant along £ x and DH does not vanish, there 
exists a number r > such that for all x £ B r (x), £ x is periodic and DH 
does not vanish along £ x . Let rj : R — > R 2 satisfy 

fl{t) = DH( V (t)) 

and 

77(0) = x. 

Then there exists 5 > such that for all t G (—5, 5) 

(a) £^( t ) is periodic and has nonzero rotation vector parallel with Q 

(b) M(j,(t))^0. 

Since -ff is constant along each orbit, (b) implies that 

6j(ti) n 6j(t a ) =0 if < *i < *2 < <5. 

and even more 

i n (h) n 4( t2 ) =0 if -5 < ti < t 2 < <5. 

See the end of the introduction section for the definition of £. 
Suppose that 

£(T) = x + (m, n) 

for (m, n) 6 Z 2 \{0}. Now consider the region Do which is bounded by £ and 
£ + (— n,m). Without loss of generality, we assume that the positions of £, 
r\ and (m, n) are as shown as in Figure 3. Let us define w in Dq as follows. 

u;|^ (t)+ (_ n>m) = r, for < t < f 
io|^ (4)+ (_ n)m ) = r - h(t), for | < t < 5 
w = 0, elsewhere in Z?o- 



Here 



2 2 
m + n 
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and h is a smooth function which is near | and r near 5. Clearly, it; is a 
smooth function satisfying 

Dw ■ V = in Do- 
Since £(s) = £(s — T) + (m, n), owing the periodicity of Cr?(t)) we have that 
as a set 

&?(t) = &?(t) + (m,n) for |t| < 5 

and 

D = D + (m,ra). 

Hence for x € Do: 

w(x + (m, n)) = x. 

Now for i e Z, we denote 

Dj = D + i(— n, m). 

Then it is clear that Dj is the region bounded by £ + i(— n, m) and £ + (i + 
1)(— n, m). Also, Dj and Dj have disjoint interior if i ^ j and 

U ieZ Di = R 2 . 

We define U : R 2 -> R as 

{/ (x) = t«(x — i(—n, m)) + rz, if x € Dj. 
Recall that r = m 2 + n 2 . Then D G (7°° (M 2 ) and 

DC/ • V = in R 2 . 
Let -S'(x) = ?7(x) — x • (— n, m). Then it is easy to see that 
S(x + (m, n)) = 5(x + (-n, m)) = S'(x). 
Hence S is rZ 2 periodic, i.e, for all »£Z 2 

5(x + rf ) = S(x). 

Now suppose p ■ Q = 0. Note that Q = ■ There exists such that 

p = u(—n,m). Let F be a viscosity solution of 

\p + DD| 2 + AV(x) ■ (p + DF) = p A . 
24 




Figure 3: Positions of £ v (t)- 

for p = v{—n,m). Assume that F(xq) — uS{xq) = mm.-g2(F — vS). Then 

\p + uDS(x )\ 2 + AV(x ) • (p + uDS(x )) > (3 A - 
Since DU ■ V = 0, we deduce that 

Pa{p) < \p + vDS{xv)\ 2 = 0(1). 

□ 



3.3 Qualitative Similarities between c p and the Asymptotic 
Growth Rate of st from RDA 



Now let us denote c*(A) as the st(A) from the RDA equation (1.1), and 
given by the variational formula (1.2). According to [39], 

c*(A) J Jn V ■ pw 2 dx 



lim 

A-5-+oo A 



sup 



J T „ w 2 dx 



where 



T = { W ^H\T n )\ V-w = 0, ||£>HI1</'(0)IHI1}- 
According to Lemma [372], it is obvious that 



Cp 2l Cp. 



More interestingly, Lemma 4.3 says that when n = 2 

lim c* = c v . 



(3.13) 



The equality (3.13) reveals that in 2d, changing the order of two limiting 
processes (one sending the flame thickness or reaction time scale to zero, 
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the other sending the turbulent intensity to infinity) does not change the 
asymptotic growth rate of the effective front speed. The following is an 
analogue of Lemma |3,4[ 



Lemma 3.5 Suppose that £ : M. — > M. is a periodic orbit with period T . 
Denote the rotation vector Q = _ Assume that Q / 0. Then 

(i) c*>0ifp-Q>0; 

(ii) c* p (A) < 0(1) forallA>Oifp-Q = 0. 

Proof: (i) Assume that p ■ Q > 0. Let x, n = n(t), 8 and (m, n) be the same 
as in the proof of Lemma 3.4 see also the Figure 3. We may choose 5 small 
enough such that for \t\ < 5 and s € M 

(£„ (t) (T + a )-£„ (t) ( S )).p>0. (3.14) 

For M > 0, we define w G PF 1,00 (T 2 ) as follows 

(w\s =M + l- 4 £ for \t\ < | 
\w = M elsewhere in T 2 . 

The definition of is at the end of the introduction section. Clearly 
Dw ■ V = 0. Also, if we choose M large enough, 

\\Dw\\ 2 <f'(0)\\w\\ 2 . 

Denote C = Ui^s^j) C T 2 , which is a closed invariant set of T 2 . For all 

Jc 
Hence 

f c V(x) ■ pw 2 {x) dx = ^ Jq f D V(£ x (s)) ■ pw 2 (£ x (s)) dxds 



c 



T fc fo V (.^( s )) • P w2 (Cx{s)) dsdx 
^ J c w 2 (x) J T ix(s) -pdsdx 



J c w 2 (x)(£, x (T) - x) -pdx 



Owing to (3.14), 



V{x) ■ pw (x) dx > M / V ■ pdx. 
c Jc 
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So according to the definition of w 

/ V{x) ■ pw 2 dx > M 2 V{x)-pdx = 0. 

Hence (i) holds. 

Now let us prove (ii). Assume p = v(—n, m) for some Due to the 



formula (1.2) 



C ;(A) < K A (p) + /'(0) (3.15) 

where k.a{p) is the unique number such that the following cell problem has 
solutions 

- Au + \p + Du\ 2 + AV{x) ■ (p + Du) = n A (p). (3.16) 

Let S be the smooth rZ 2 periodic solution from the proof of Theorem 3.4 
with r = m 2 + n 2 . Assume that 

u{xq) — uS(xo) = min(-u — vS). 

Combining with the equality that V(x) ■ (p + uDS) = 0, we have that 

ka(p) < -vAS(x ) + \p + vDS{x Q )\ 2 = O(l). (3.17) 

□ 

Lemma 3.6 Let £ : R — > M 2 be a periodic orbit with rotation vector Q. 
Then there must exist another periodic orbit £ with rotation vector XQ for 
some A < 0. 

Proof: This is obvious if Q = 0. So let us assume that Q / 0. Since n = 2, 
any two periodic orbits must have parallel rotation vectors. So according 
to the Poincare recurrence theorem, for a.e x, £ x is a periodic orbit with 
rotation vector parallel to Q. Hence for a.e x, 

ton Z*F)- X =Xx Q. 

for some \ x £ M. Let £(0) = x. Then = 1 and when x is close to x, ^ x is 
also periodic and X x > 0. Note that 

= J J2 V(x) -Qdx =±Jq J T2 V(£ x (t)) ■ Q dxdt 

= T f T 2foV(£ x (t))-Qdtdx 
_ r (S*(T)-x)-Q r 
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Sending T —> +00, by dominant convergence theorem, we derive that 

X x dx = 0. 



IJ2 

So there must exist x such that is periodic and A^ < 0. □ 

The following theorem says that the G-equation model, the F-equation 
model and the reaction-diffusion equation model predict the bending effect 
simultaneously, i.e, if one of these three models predicts the bending effect, 



so will the other two. It follows immediately from Lemma 3.4, Lemma 3.5 



(3.17) and Lemma 3.6 



Theorem 3.3 Assume n = 2. Then one of the following holds, 
ft) 

Cp = c* p = 0, for all peM 2 . 
(ii) There exists a nonzero vector Q S M 2 such that 

c p > c* > 0, for all p such that p • Q 7^ 

and for all A > 

aA{p),f3A(p),7A(p),K A (p),c* p (A) < O(l), if\p\=l andp-Q = 0. 

Moreover, case (ii) happens if and only if there exists a periodic orbit with 
rotation vector Q / 0. 

For reader's convenience, we recall the notations. 

ctA(p) the turbulent flame speed predicted by the G-equation 
1a{p) the turbulent flame speed predicted by the F-equation 
^c p (A) the turbulent flame speed prediced by the scalar RDA equation 

and 



j (3a(p) effective Hamiltonian from the F-equation (3.2) 

1 k a(p) effective Hamiltonian from the viscous F-equation (3.16). 



Also 



°A(p) _ . Pa(p) _ 7A(p) 



lim — - — = lim — - — = lim — - — (Theorem 3.1) 
A^+oo A A A^f+oo A ' — ' 
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and 



hm — - — . 

A^+oo A 



A good example to demonstrate the above theorem is the following cat's 
eye flow. 

Example 3.1 (cat's eye flow) For the cat's eye flow, the stream function 
is H = sin2-7ra;i sin27TX2 + 5cos2-7ra;i cos2-7rx2 for 5 6 [0, 1]. Note that when 
5 = 0, it is the cellular flow which belongs to case (i) in Theorem 3.3. For 
5 > 0, the zero level curve {H = 0} is a periodic orbit and has a rotation 
vector parallel to (1,1). So it is the case (ii) in Theorem\3.3\ Hence 



c p , c* > 0, if p is not parallel to 



-1,1) 



and for all A > 

ola{p), Pa(p), 7a(p)> ka(p), c* p (A) < 0(1), if p is parallel to (-1,1). 

Hence for 5^0, the bending effect only occurs along the direction parallel 
to (-1,1). 



3.4 Growth Laws of St from Different Models in Cellular 
Flows 



In the case (i) of Theorem 3.3 it remains an important and challenging 
problem to determine the exact growth law of ctA, 1A and c*(A) as A — > +oo. 
When V is the cellular flow, it is known that 



a A = 



A 



log A 



ifp=(l,0) 



and 



c*JA) = 0(Ai) ifp=(l,0). 



(3.18) 



(3.19) 



Due to (3.3), (3.4), (4.4) and (1.2), it is easy to see that cca(p), Pa(p), 



1 A, k a(p) and c*(A) are all convex as functions of p. The symmetry of 
the stream function H also implies that these four functions are even and 
symmetric: f(p) = f(—p) and 

f((Pl,P2)) = f((P2,Pl)) = f{{-Pl,Pl))- 

Note that for ka and c*(p), we also need to use the fact that two adjoint 
compact operators have the same eigenvalues. Hence (3.18) and (3.19) are 
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true for all p 0. In the following, we show that Pa and 7A also enjoy the 
y^gAj law. The exact growth law of ka for cellular flow remains an open 
problem. See Remark 4.1 for more details. 

Theorem 3.4 IfV is the cellular flow, then 

ia(p),Pa(p) = o' 1 

for all \p\ = 1. 



\og(A) 



Proof: Throughout this proof, C represents a constant which only depends 
on V. According to Lemma |3.1| and the above discussions, it suffices to 
establish this for (3 A (p) and p = (1,0). It is not hard to show that for any 
two points x, x £ Qi = [0,1] x [0,1], there exists a Lipschitz continuous 
curve f : [0, L) R 2 such that f (0) = x, f (L) = x and 



|f(t) - VAV(£(t))\ < 1 for a.e t G [0,L]. 
Moreover L < C. Let U(x) = p- x+ ^ x ~> for A > 1. Then {/ satisfies that 

\DU\ 2 + s/AV ■ DU = ^ < C. 

A ~ 

Hence 

\U(x 1 )-U(x Q )\<C. 
So if we choose F such that J T2 F dx = 0, then 

|F| < C(vCi+l). 



Hence due to Step 2 in the proof of Lemma 3.3, we obtain that 

\p + DF\ < C{A^ + 1). (3.20) 

Write \p + DF\ and denote c\a as the effective Hamiltonian 

of the following modified G-equation 

co A \p + DG\ + AV{x) ■ (p + DG) = a A . 
Now we claim that 

Pa < aA- 
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In fact, suppose that <p £ C 1 (T n ) and 

(p(xo) — F(xq) = min(0 — F). 

Then 

\p + D(f>(x )\ 2 + AV(x ) • (p + Xtyfro)) > /3a- 
Also, it is easy to see that 

|p + 4>(x )\ < UJ A - 

Hence 

w A |p + D(j)(xo)\ + AV(x ) ■ (p + Dcf)(x )) > p A . 

So 

max{u; A \p + D0| + AV ■ (p + > /3 A - 

Owing to the inf-max formula ( |3.3[ ), my claim holds. Hence 

A /A 

PA < OLA = W A a^ < C- -tt : r < O -r-rr 

-a log(A) - log(uu) V lo g(A) 



The last inequality is due to (3.20). 



Summary: The growth laws of the turbulent flame speeds from different 
models in the cellular flow are as follows. For p ^ 0, 



<*A(p) = O yj^Ay 
7A(P) = O (i^) 
{c* p (A) = 0(A\). 



A remaining interesting question is what is \im A ^ t+00 (7a~ ola) or lim A _j. +00 ^ 
for the cellular flow. We will investigate this in a future paper. 



4 Viscous G-equation and Viscous F-equation 

Recall that without loss of generality, we set d = r r = si = 1. For A > 0, we 
denote \A as the effective Hamiltonian of the cell problem associated with 
the viscous G-equation 

-AG+\p + DG\ + AV(x) ■ (P + DG) = X A, (4.1) 
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and ka as the effective Hamiltonian of the following cell problem 

- AS + \p + DS\ 2 + AV(x) • (p + L>S) = ka- (4.2) 

Here we omit the dependence of xa and ka on p. As it was mentioned in the 
introduction part, \A can be viewed as a simplified model of the turbulent 
flame speed predicted by the curvature G-equation (1.10). Also, ka is a 
simple upper bound for the turbulent flame speed st predicted by the RDA 
model and the Majda-Souganidis model for a = 1. Precisely speaking, owing 
to (jTJ, 

c* p (A) < « A + /'(0). 

The main purpose of this section is to identify the asymptotic growth rate 
of ka in two dimensions. Similar to the inviscid case, both xa and ka are 
given by inf-max formulas: 



XA= inf max{-A</> + \p + D(j)\ + AV(x) ■ (p + D6)}, (4.3) 

0eC 2 (T«) T n 



and 



«A = inf max{- A</) + \p + Dcj)\ 2 + AV{x) ■ (p + D<f>)}. (4.4) 

0eC 2 (T n ) T n 

Clearly, xa, ^A < O(A) and xa < K A + j- The following theorem simply 
says that the asymptotic behavior in the viscous case can not be larger than 
that in the inviscid case. 

Lemma 4.1 

limsup — < limsup — < c p . 

A— >+oo A A— 5>+oo A 

Proof: The first inequality is obvious. Let us prove the second one. For 
e > 0, by routine mollification, there exists a <p £ C 2 (T n ) such that 

maxjn{V(x) ■ (p + D(p e )} < inf^ eC i( T n) maxjn{V(x) ■ (p + Dcf))} + e 

= c p + e. 



Then by (4.4), 



— f- < — max(-A0 e + \p + D<p e \ 2 ) + c p + e. 
A A T n ^ 



So the inequality holds by sending A — > +oo. □ 
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Due to the weak nonlinearity in the viscous G-equation, lim J 4_> +00 xa/A 
might be smaller than c p . See [TH] for the example of shear flows. 

The following theorem says the equality holds for ka in two space di- 
mensions. Precisely speaking, 



Theorem 4.1 For n = 2, 



lim — - = Co. 



We first prove several Lemmas. 
Lemma 4.2 



lim inf — ^ > sup / p ■ V da, 



where A is the collection of smooth invariant measures, i.e, 

A = {a = w 2 dx\ w G fl^T 71 ), / w 2 dx = 1 and V(x) ■ Dw = 0}. 



Proof: Choose w G A. Multiply tu 2 on both sides of equation (4.2) and 
integration by parts. We deduce that 



— / wDSDw dx + — [ \p + DS\ 2 w 2 dx + / p-V(x)w 2 dx 

A Jjn A Jjn Jfn 

By Cauchy inequality 

| f wDSDwdx\<-([ \DS\ 2 w 2 dx+ [ \Dw\ 2 dx). 

JT™ ^ Jfn Jfn 

Sending A — > +oo, we have that 

lim inf — ^ > / p ■ V(x)w 2 dx. 

A-^+OO A Jfn 

□ 



Recall that Lemma 3.2 says that 



Cp = max / p ■ V(x) da, 

crEA Jfn 



where A is the collection of all Borel proba bilit y m easu re on T n invariant 
under the flow £ = V{^). Combining Lemma 
immediate corollary. 



4.1 



and 



4.2 



the following is an 
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Corollary 4.1 // 



sup / p ■ V(x) da = max / p ■ V(x) da. 



(4.5) 



Then lim^_ 



K-A 

-oo a 



V ■ 



As we mentioned after (3.13) in section 3.2, the physical meaning of 



(4.5) is that changing the order of limits of sending the flame thickness (or 



reaction time scale) to zero and sending the turbulent intensity to infinity 
does not change the front speed asymptotic growth rate. It is not clear to us 



when (4.5) holds. Nevertheless we are able to establish it in two dimensions. 



Lemma 4.3 When n 



sup / p ■ V(x) da = max / p ■ V(x) da. (4-6) 



Proof: If c p = 0, then (4.6) is obvious. So let us assume that c p > 0. Then 



due to Corollary 3.1 there exists a periodic orbit £z : [0,T] —> M? such that 
£2(0) = x and 



p-(&(T)-x) 







p ■ V(x) da 



where T is the period and a the Borel probability invariant measure defined 
as 

/ /(*) da = ~ [ f(m) dt for all / G C(T 2 ). 
Ji 2 1 Jo 

The strategy of the proof is to construct a sequence of smooth invariant 



measures to approximate a which implies (4.6). Let 5, rj = rj(t) be from the 



proof of Lemma 3.4 See also Figure 3. Choose a sequence of functions {4>n\ 

such that for all n G N, 4> n G C°°(T 2 ) and 

(i) 



(ii) 
(i") 



'&j(t) 



ct for some q G M; 



in T 2 \U w< i4 ( t); 
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(iv) 

Then it is clear that 



„ dx = 1. 

T 2 ' " 



V[x) ■ D<t> n = 0. 



Let us assume that 

j.2 



a weekly as Random measures on T™. 
Then a must an invariant probability measure supported on £. Hence 

a = a. 

□ 



Proof of Theorem \4.1\ It follows immediately from Corollary 4.1 and 
Lemma 14.31 



Remark 4.1 It remains an open problem whether lim^^+oo ^ exists for 
n > 3. We also do not know the exact growth law of ka in cellular flow 
except for the sublinearity. Since c*(A) < ka + /'(0), ka has a growth lower 

bound of 0{A~i). However, it is proved in \2(tf that xa is uniformly bounded 
as A — > +oo, i.e, for all A > 

^M<o(i). 
\p\ 

A remaining challenging and very interesting question is whether the origi- 



nal curvature equation (1.10) will predict similar uniform boundness of the 
turbulent flame speed. The physical meaning is whether the curvature effect 
could have essential contributions to the "strong bending effect" observed 
in experiments |31j . 



5 Appendix: Review of Viscosity Solutions of 
Hamilton- Jacob i Equation and Cell Problem. 

Suppose that is an open set in W 1 and H = H(p, x) G C(IR n x Q). Let us 
look at the following Hamilton-Jacobi equation 

H(Du,x) = inn. (5.1) 
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It is well known that the above equation might not have classical (C 1 ) solu- 
tions. Solutions need to be interpreted in the viscosity sense. See the User's 
guide [10J for reference. 

• Viscosity subsolution (supersolution): We say that u G C(ft) is a 
viscosity subsolution (supersolution) of 

H(Du,x) = in ft 

if for <f) G C 1 (ft) and xq G ft satisfying that 

4>(xq) — u(xo) < (>)4>(x) — u(x) for x G ft 

then 

H(4>(x ),x ) < (>)0. 

• Viscosity solution: We say that u G C(ft) is a viscosity solution if it is 
both a viscosity subsolution and supersolution. 

• Approximate viscosity solution: We say that the equation 

H(Du,x) = 

admits approximate viscosity solutions if for any e > 0, there exists u e G 
C(ft) such that it is both a viscosity subsolution of 

H(Du e ,x) = e 

and a viscosity supersolution of 

H(Du e ,x) = — e. 



• Connection between viscosity solutions and classical solutions: 

Any C 1 solution of equation (5.1) must be a viscosity solution. If u is a 
viscosity solution and is differentiable at x G ft, then it satisfies the equation 
at x, i.e, 

H(Du(x),x) = 0. 



• Coercivity of H and Lipschitz continuity of viscosity solutions: 

Suppose that H is coercive in the p variable, i.e 

lim H (p, x) = +oo uniformly in ft, 

\p\— >+oo 
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then any viscosity solution of equation (5.1) is W 



rl.oo 



• Cell problem, Homogenization and the effective Hamiltonian: 

Suppose that H 6 C(IR n x T n ) and is coercive in the p variable. It was 
proved in [18] that for any P £ W 1 , there exists a UNIQUE number H(P) 
such that the following equation 

H(P + Du,x) = H(P) inIR n 

admits periodic viscosity solutions. This is the so called "cell problem" in the 
homogenization theory. The quantity H{P) is called effective Hamiltonian. 
For e > 0, assume that u e is the unique viscosity solution of 

[u\ + H(Du e , |) = 
\u e (x,0) = g(x). 

Then as e — > 0, u e uniformly converge to the u which is the unique solution 
of the following effective equation 

(u t + H(Du) = 
|u(x,0) = g(x). 

The cell problem and homogenizaton analysis were performed for fully non- 
linear PDEs in f!2l. 



• Inf-max formula for convex effective Hamiltonian: If H is also 
convex in the p variable, the H{P) is given by the following inf-max formula. 
See [15J for instance 

H(P)= inf max if (P + DS.x). 

0ec 1 (ir n ) T" 



6 Concluding Remarks 

We analyzed the turbulent flame speeds (st) of G-equation, the analogous 
quadratically nonlinear Hamilton-Jacobi equation (the F-equation) and the 
reaction-diffusion-advection equation in the presence of steady and periodic 
compressible flows, shear flows, and incompressible flows (cellular and cat's 
eye flows). The strain effect slows down flame propagation in Id compress- 
ible flow and shear flows. If f'(0) is no less than a threshold value, the G- 
equation always predicts smaller turbulent flame speeds than the F-equation 
model. The F and G equations share the same asymptotic growth rates in 
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incompressible flows, in particular, the sy's obey the same sublinear growth 
laws i 4 n in cellular flows. Moreover, in two dimensions, the three mod- 

log(A) ' ' 

els predict the bending effect simultaneously. Moreover, in 2d, changing 
the order of limits of sending the flame thickness (or reaction time scale) 
to zero and sending the turbulent intensity to infinity does not change the 
asymptotic growth rate. The nonoccurrence of the bending effect is charac- 
terized precisely by the existence of periodic orbits on the torus and their 
rotation vectors in the underlying advective flow field. The effective Hamil- 
tonian of the viscous F-equation has the same growth rate as the inviscid 
F and G-equations in two dimensions. Asymptotic growth rate of the effec- 
tive Hamiltonian from the viscous F-equation in higher dimensions will be 
explored further in a future work. 
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